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Abstract 



A particular class of variant axion models with two higgs doublets and a 
singlet is studied. In these models the axion couples either to the u-quark 
or t-quark or both, but not to b, c, s, or d. When the axion couples to 
only one quark the models possess the desirable feature of having no domain 
wall problem, which makes them viable candidates for a cosmological axion 
string scenario. We calculate the axion couplings to leptons, photons and 
nucleons, and the astrophysical constraints on the axion decay constant Va 
are investigated and compared to the DFSZ axion model. We find that the 
most restrictive lower bound on Va, that from SN1987a, is lowered by up 
to a factor of about 35, depending on the model and also the ratio of the 
vacuum expectation values of the higgs doublets. For scenarios with axionic 
strings, the allowed window for Va in the u quark model can be more than 
two orders of magnitude. For inflationary scenarios, the cosmological upper 
bound on Va/N, where N is the QCD anomaly factor, is unaffected: however, 
the variant models have N either 3 or 6 times smaller than the DFSZ model. 
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I. INTRODUCTION 



The relevance of instantons |T[| to physics was first shown by 't Hooft [^] in his resolution 
of the f/A(l) problem It was first pointed out by him that the topological structure of 
the vacuum of any non-abelian gauge theory, which includes QCD, is non-trivial. There 
are gauge transformations, characterized by different topological numbers n (the Pontryagin 
index), which can not be continuously deformed to one another. This gives rise to distinct 
ground states, labelled by different n and separated by finite energy barriers. Instantons can 
be physically interpreted as quantum mechanical tunneling events, in Euclidean spacetime, 
between these different ground states One then has to construct the true, gauge- 

invariant, vacuum out of these degenerate n- vacua by taking a linear combination 

\e)=Y.e-'-%n). (1) 

n 

In the path integral approach, the effect is to add the so-called 6'-term, to the ordinary QCD 
lagrangian, or 

Under the combined action of charge conjugation and parity transformation the 6'-term 
changes sign and hence it violates CP invariance. There is also CP violation communicated 
from the quark mass matrix M: if we diagonalise it with a bi-unitary transformation, we 
find that the coupling constant 9 is modified to 

^ = ^ + arg(detM). (3) 

Despite being the coupling constant of a total derivative term, the parameter 6 is observable, 
through its effect on the neutron electric dipole moment . The current experimental upper 
limit 1^ on the electric dipole moment implies that 6 is less than about 10~^. 

Such a small value of 6 contradicts our expectation that a dimensionless free parameter 
should be of order one. Many ideas have emerged in trying to resolve the puzzle. One 
of the most elegant solutions was proposed by Peccei and Quinn in 1977 The authors 
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postulated the invariance of the lagrangian under the transformations of a new extra global 
chiral U (1) symmetry, called PQ-symmetry, thus enlarging the symmetry group of the SM to 
SU (3)c X SU (2) X f/(l)y X f/(l)pQ . To accomodate the extra charges of the new symmetry 
one needs (at least) one extra higgs doublet. When this symmetry is spontaneously broken, 
a pseudoscalar boson appears ^ in the theory, called the axion. Normally, one would think 
that the axion is massless, as it is the Nambu-Goldstone boson of the PQ symetry. However, 
the U{l)pQ symmetry is an anomalous one, spoiled by the effect of instantons in the QCD 
vacuum, and the axion picks up a small mass via the axion-gluon-gluon triangle anomaly. 

The assignment of appropriate PQ-charges to the higgs fields and consequently to the 
quarks is responsible for the presence of the anomaly in the PQ current, and also for the 
variety of the axion models. In the original model, known as the PQWW (Peccei-Quinn- 
Weinberg-Wilczek) axion model 0, all the quarks of the same chirality were assigned the 
same PQ charge. Unfortunately the model was ruled out both by particle physics exper- 
iments and by astrophysical observations. The former gave constraints which came from 
J/ip and T meson decays, reactor and beam dump experiments and nuclear deexcita- 
tions [0. The latter are in fact more restrictive and imply that Va > 10^° GeV for most 



axion models considered to date |TT|. A way out of this was first proposed by Kim in 1979 



and subsequently by Shifman, Vainstein & Zakharov ||12|, who added a higgs singlet, thus 
enabling the axion decay constant to be much higher than the electroweak scale. However, 
it cannot be too high, as it is also possible to restrict Va from above through cosmological 
arguments. Coherent oscillations in axion field can be produced after inflation or via the 
formation of axionic cosmic strings. The requirement that the energy density in these oscil- 
lations is not large enough to overdose the universe puts an upper limit on Va- The current 
values on these limits for the axion decay constant are Va less than about lO^^GeV |T3|JTl 



from inflationary scenarios and Va ^ 2.6 x lO^^GeV [0] from axion strings (for Hq = 50 
km s~^ Mpc~^, and making the conservative assumption that radiation from infinite strings 
dominates that from loops). As we can see there is only a very small window left for the 



axion. 



Cosmology also restricts on the value of A^, the parameter that characterises the QCD 
anomaly, which is related to the number of quarks that couple to the axion. If there is 
no inflation between the PQ symmetry-breaking transition and the present day, a dense 
network of axion strings is formed. At around a temperature of 1 GeV, each string becomes 
the junction of (in our normalisation convention) 2N domain walls [|r^. In order to avoid 
the domain walls dominating the energy density, 2N must be equal to unity, so that the 
string-wall system can annihilate. Models with N > ^ must have a period of inflation at 
a low energy scale, or must reheat after inflation to less than the PQ symmetry-breaking 
temperature, to remain viable. 

In this paper we examine some variant axion models based on those proposed first by 



Peccei, Wu & Yanagida [|T6| and independently by Krauss & Wilczek [0. Their models 
were constructed with two higgs doublets and assigned of different PQ charges to different 
quarks. The original reason for this is that in order to make an axion model which avoided 
the particle physics constraints at the time it was essential to decrease the axion couplings 
to c and b quarks on one hand, and on the other to sufficiently suppress the tt+o 
decay rate. One must also have a sufficiently short lived axion so it cannot be detected by 
the other experiments. To accomplish that, one has to couple both the c and b quarks to 
the same higgs field. However, as we know the limits on the strangeness changing neutral 
currents are very tight and for that reason we have to couple the strange quark to the same 
higgs doublet as c and b. Thus the only way to realise the Peccei-Quinn symmetry is to 
couple either the u or the t or both quarks to a second higgs doublet. So we have three 
different models in our hands, two of which have the cosmologically desirable property that 
N = ^, and hence have no domain wall problem. 

Although the original models are also ruled out, along with the PQWW axion, by the 
astrophysical constraints, they also have extensions with a higgs singlet. The astrophys- 
ical constraints on these models, presented in the next section, have not previously been 
considered, and we present an analysis in this paper of the bounds on the axion scale that 
can be inferred from their couplings to electrons, photons, and nucleons in astrophysical 



processes. We find, as for the standard DFSZ |T8[ and KSVZ [jT2[ invisible axions, tfiat 
tfie tigfitest constraint comes, via tfie nucleon coupling, from SN1987a. However, in these 
variant models, the coupling is generally weaker, and weakens the lower bound on Va by a 
factor of between 1.4 and 35, depending on which quarks the axion couples to, and on the 
ratio of the vacuum expectation values of the higgs doublets. The cosmological upper bound 
on Va from the axion density is also reduced, by a factor of either 3 or 6, as a result of the 
smaller QCD anomaly factor. 

II. DESCRIPTION OF THE MODEL 



The models we are going to discuss were first proposed by Geng & Ng [jl9[ and has 
elements from both the previously discussed one and that of DFSZ. There, an extra higgs 
singlet (f) is introduced, the phase of which is reserved for the axion. A direct coupling of 
(j) to quarks and leptons is impossible; but in the DFSZ model it couples to both (pi and 
02- We however couple only to one of the two doublets, namely 0i, which then couples 
to the 'special' quarks according to the three models discussed above. The other higgs field 
couples to the rest of the quarks and the leptons. Therefore we have the following Yukawa 
interactions: 



C-Y = mu<pium) + fi{qu(p2dRj) + ^ frjilu4>2UR,) + h.c. (Model I) 

i=2,3 

J^Y = f^siqiihum) + ftjiqii^dRj) + ^ fijiquhuRj) + h.c. (Model II) 

i=i,2 

'^Y = fii{qLi(t)iURi) + fiz{qLi(t)iUR^) + ffjiquhdRj) + fi2{qLihuR2) + h.c. (Model III) 



where (j) = ia24>*. Our nomenclature is the same as that of Peccei, Wu & Yanagida [IB 
In the first model it is the u quark that couples to the axion, whereas in the second is the 
t. Finally in the last model both quarks couple to the axion. 

The most general renormalizable potential for the model, consistent with gauge, as well 
as PQ symmetry and renormalisability is [|l^ ,|T9|l 

V = Ai(0l0i - vl/2y + X2{4h - vl/2y + A(0*0 - t;V2)2 + 



(a</)l(/)i + b(j)l(j)2)(f)*(j) + c{(f)^icr2(f)24>'^ + h.c) + d\(j)^ia2(j)2\'^ + e\(f)\(f)2^^ 



(4) 



The appropriate PQ transformations according to which the quarks acquire a PQ charge, 
as well as leaving the Yukawa lagrangians invariant are 



iQflja 



oi{QLj-QRj)oi 



(5) 



where a is explained below. In the models under discussion, we can choose an assignment 
of PQ charges such that Qlj = and only some of the Q^j ^ 0. Furthermore, we impose 
the normalisation condition that Q^j = 1 for the 'special' quarks, which are ur and/or t/j. 
That further fixes the transformations for the higgs fields, which are as follows: 

1 



= -j=vi exp [-i{a + az)], 
[h) = -^V2exp (iaz), 
{(p) = -^vexp (ia), 



(6) 



where Vi, V2, v are the vacuum expectation values of the higgs fields and a, az are the angles 
conjugate to the axion and the longitudinal degree of freedom of the boson respectively. 

Our following step is to find an expression for the axion decay constant. Let a' and Z 
be the two goldstone bosons after the breaking of the Peccei-Quinn symmetry and before 
instanton effects are taken into consideration. The first one, a', is the massless axion and 
the second, Z, the goldstone boson that is eventually eaten by the Z^. One then has the 
following equation 



( 



V 







( \ 



a 



(7) 



The 2x2 matrix on the right hand side of (|^) is the most general matrix compatible with 
the requirement not to mix the axion with the Z^ boson. On the other hand if one takes 
the kinetic term for these massless scalar degrees of freedom then one has 
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(8) 



The derivative of (|^) compared with (j^) gives the expression for the axion decay constant 
as well as the electroweak breaking scale 



Va 



\ 



VIV2 

2 i 2 ' 



vew = vii = Jvf + f i = 246 GeV, 



fio 



(9) 
(10) 
(11) 



The expression is the same for all three models and as we can see, it is v that ultimately 
fixes the scale for Va- 

On the other hand if we apply the same normalisation convention to the DFSZ axion, 
that is if we assign Q/jj = 1 to every quark, then the transformation of the higgs fields are: 



bi) = -j=vi exp [-i{a + az)], 
62) = -^V2 exp [i{-a + az)], 



(12) 



1 

72 



V exp [la) 



where 



2 \ 2 ' 



vew = Vu 



2 - 246 GeV, 



(13) 

(14) 
(15) 



Our choice of charges ensure that f^, which is essentially equal to the PQ symmetry-breaking 
scale, is also equal to the axion decay constant fa (as defined by Srednicki ||2l|]). We have 
avoided using because of its many definitions in the literature. 
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III. THE AXION COUPLINGS 



Our next step is to determine the couplings of the axion to different quarks and leptons. 
In the basis where quarks q'^ are eigenstates of the weak interaction, the QCD part of the 
lagrangian is 

^QCD = qjiilf,D^)qj - qLiMijq'nj - qmMjjqLj - -G'^'^Gb^^ + 9^^G'^'"Gb^,u, (16) 

where M is the quark mass matrix, and depends on a and az, and hence a'. If we diagonahse 
the quark mass matrix, we find that in the mass eigenstate basis the lagrangian is 

JC-qcd = qA^l.D^ - m^)q, - ^G^'G^,, + ^{d.a'f + " 2N^JG^,'Gb,. - 

-^fa-q,^,^,q,, {11) 

where N is explained below. If there were no mixing between the axion and the the 
coupling Qj would just be the difference of the PQ charges of the right and left-handed 
quarks qj. The presence of the mixing modifies the relation to 

= {Qr, - Ql,) - {Ynj - Yi^.f-^, (18) 

where Y^j and Y^j are the hypercharges of the right and left chiral fields respectively of the 
j-th quark. 

The constant depends on the number and type of particles with the gluon anomaly 

N = ^{Qnj-QLjh, (19) 
j 

where tj is the index of the SU{3)c representation to which the fields belong (for the known 
quarks t = |). Hence, in our normalisation convention, for Models I and II, = i, and for 
Model III, = 1. As advertised, the first two models have no axionic domain wall problem. 
In the case of the DFSZ axion we have = 3. 

When calculating the axion-quark interactions in practice, the derivative interaction in 
(P!7|) can prove troublesome, and it is usual to leave the phases corresponding to the axion 
degree of freedom in the quark mass matrix, so that the interaction term is 
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= -m,-^,e^^^(^'2.-+Ay,'..)^.^ (20) 
where AQj — Qrj — Qlj and AYj — Y^j — Ylj. To first order in a' /va-, tiiis gives 

Tfl ' 

l^aq = -Wj — (l'qjj5qj- (21) 

The (pseudoscalar) Yukawa coupling of the interaction between the axion and the j-th quark 
is then clearly 

Kj = 9j-f. (22) 

The interactions with leptons can be calculated similarly. The term involving the axion and 
the Z-th lepton ei is 

>Cai = -igi—a'ea^eu (23) 
from which we define the Yukawa coupling 

hai = 9i—. (24) 
It is useful for later calculations to express the couplings in terms of tan /3, where 

tan/3=— . (25) 



So for Model I: 



tan^/S 



where o?j — d,s,b and = e, r. Using the same notation we can work out the couplings 
for Model II, 
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and for Model III: 



gu,t 



9dj ,ei 



9c 



vl 



vf + vl 1 + tan^ (3 
1 



sin^ P, 



vf + vi 1 + tan^ p 
1 



cos^ P, 



cos^ p. 



(32) 
(33) 
(34) 



vf + vi 1 + tan^ p 
Lastly, we write down the interaction between the axion and photons, which arises in 
the same way as for the axion-gluon interaction, through the 077 triangle diagram 



where 



r = N a'F F^" 



(35) 



(36) 



In this equation, Cj is the electric charge of the j-th fermion (with each colour quark counted 
separately) . 

It is important to note that below the QCD scale (~ 200MeV), free quarks do not exist 
so one has to consider the effective couplings of axions to nucleons, which arise from the 
axion mixing with 7r° and the rj' . It is exactly this mixing that allows us to calculate the 
axion mass and also has an effect on the coupling to two photons. 

Standard current algebra techniques p0| , |2l[| tell us that the physical axion a has mass 
(including the strange quark) 

1/2 



rria = N 

Va 



and that its coupling to photons is 



'l + z){l + z + w) 



(37) 



where 



47r2 



r = -h aFF 



Va [n 3(1 + Z + w) 



(38) 



(39) 
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and we define z and ti? as 



Tfl 

z= — = 0.553 ± 0.043, 



Tfl 

w = — = 0.029 ± 0.0043. 



For models I & II one has {N, N^) = (1/2, 4/3) and for model III (A^, N^) = (1, 8/3) . As we 
can see the ratio Ng/N = 8/3 for all of them. This ratio holds both for the DFSZ and the 
KSVZ axion, as well as for GUT axion models based on SU(5) ||21|| . 

The same current algebra techniques enable us to compute the axion-nucleon coupling 
|2T|. First we write down the anomaly-free axion current, including the strange quark 

1 / 2 , \ 1 / 2z 



3 a = Vadf,a + -(gu- — ■ — n] u-if,-i^u + -(gd- -rV~^ — ^) ^7^75^^ + 

^ 2V l+z+w J 2V 1+z+w J 

+ 9 - 1 ^ ^ ^ ^7m75S. 40 
2\\+z+wJ 



To find the couplings with the nucleons, we take the axionic current in the effective theory 
of nucleons, 

= ^^(^0 + ^3^3)7^75^, (41) 
where is the usual Pauli matrix, and ip is the nucleon doublet 



sandwich it between nucleon states, and compare with the expression obtained by using Eq. 
(PH). We find that the isoscalar and isovector couplings are given by 

go = ]:{Au + Ad- 2 As) (g^ + ga - 2g, - i±^^^2Ar) + 
6 V I + z + w J 

+-{Au + Ad + As){gu + gd + gs-'^N), (42) 
^73 = i {Au - Ad) [g^ -gd- ^|~^^ 2Ar) , (43) 

where s^Ag = ('^|^7^755'|V'), with being the spin of the nucleon. We can find values for 
the Ag's by combining the measurement of the proton spin structure function by the E143 
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collaboration p3[, which gives Am + Ad + As = 0.27 ± 0.10, with quark model and SU(3) 



flavour sjTiimetry predictions for nucleon and hyperon couplings. We flnd 

Au = 0.813 ± 0.072, 
Ad = -0.444 ± 0.072, 
As = -0.10 ±0.04. 

Translating these results into the language of the effective lagrangian approach we find 

Tn I — 

^^ai>^ = -i «V'(^o + ^3^3)75^- (44) 

It is now straightforward to find the Yukawa couplings for the axion-proton and axion- 
neutron interactions. For Model I, 

TTl 

hap = —[-(0.678 ± 0.038) cos 2/3 - (0.226 ± 0.058)], (45) 

in 

han = —[(0.578 ± 0.038) cos 2/3 + (0.136 ± 0.058)], (46) 



for Model II 



in 

hap = ^[-(0.678 ± 0.038) cos 2/3 - (1.039 ± 0.058)], (47) 
h^^ = ^[(0.578 ± 0.038) cos 2/? + (0.579 ± 0.058)], (48) 



and Model III 



Tfl 

hap = —[-(0.678 ±0.038) COS 2/3- (0.587 ±0.096)], (49) 

TTl 

h^^ = — [(0.578 ±0.038) cos 2/? ± (0.137 ±0.096)]. (50) 

Va 

Finally the Yukawa couplings for the DFSZ axion are 

TTl 

hap = ^[-(1.357 ±0.077) COS 2/3 - (1.896 ±0.275)], (51) 

Tfl 

han = — [(1.157 ±0.077) COS 2/3 ± (0.276 ±0.275)]. (52) 

"^a 
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IV. ASTROPHYSICAL CONSTRAINTS 



We are now interested in seeing how the astrophysical hmits serve to constrain the 
couphngs and hence the PQ scale Va- In our case these constraints come exclusively from 
the application of energy loss arguments to stars ||lTl. According to it, if there are low mass 
particles like neutrinos or novel particles interacting weakly with matter and radiation like 
axions, they can be produced in large numbers in stellar interiors and can afterwards escape 
freely. In this way the stars are drained of energy and so alter their standard evolutionary 
course. 

We can now proceed to the following step, which is checking how these three models 
behave under constraints coming from astrophysics. It would also be useful to compare 
the results with the DFSZ model, mainly for two reasons. First of all, it has significant 
similarities with our models, especially in the way that fermions acquire their PQ charges. 
Secondly, it is a well studied and established axion model and all existing constraints refer 
mostly to it. 

We begin by estimating the bounds on the axion couplings to electrons. The most 
restrictive bound comes from the so called helium ignition argument m red giants |1TT| , |24| , 
which states that if a red giant produces a large number of neutrinos or other weakly 
interacting particles, then helium will ignite in a much later time because larger density 
will be required. For sufficiently light particles the energy loss rate would be so large, that 
helium would not ignite at all, so the stars would directly become white dwarfs. The limit 



is HI 



hae < 2.5 X 10~". (53) 

So according to equations (p^, (pT]), (pO|), (|33| ) and ([53| ) we get the following lower limit for 
the axion decay constant: 

Va>2yi 10^ cos^ (3 GeV, (54) 

which holds for all models. Another important interaction axions have is the one with 
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photons. For the models under discussion it was first examined by Cheng, Ceng fc Ni p5 . 
The hmit set on the constraint for the axion-photon coupling is ^ 



ha^<lx 10"". (55) 
From this and equation (24) is straightforward to conclude that 

Va > 8.5 X 10^ GeV, (Models I, II) (56) 

and 

Va > 1.7 X 10^ GeV, (Model III) (57) 

The strictest bound though arises from the axion coupling to nucleons. The most restrictive 



value for this constraint come from the supernova SN1987a neutrino signal [|2^. The general 
argument underlying these kinds of measurements is as follows. We assume that we have 
some new light particles, produced in the interior of a neutron star, more weakly interacting 
than neutrinos e.g axions, and their Yukawa coupling noted as ha^, where ip is the nucleon 
doublet. There are two possibilities. If ha^p is too small, then axions cannot be trapped in 
the interior of the star. Their mean free path is bigger than the radius of the star and so 
they can escape freely. In this case axions are produced inside the whole volume of the star 
and so the axion flux would be La oc /i^^. The reason for this is that axion production would 
be dominated by processes like axion bremsstrahlung from nucleons i.e n + p — > n + p + a. 
Hence, La increases with increasing ha^p and eventually the axion flux equals the neutrino 
flux L^. Let this value of the coupling be hmin- On the other hand, for larger /i^^, then 
axions can be trapped and thermalised, in which case they are emitted from a sphere of 
radius Ra- The stronger the coupling, the larger Ra is. With blackbody surface emission. 
La oc R^T'^^Ra), where T{Ra) is the temperature at radius Ra- For a nascent neutron star, 
R'^T'^{R) is a rapidly decreasing function of R. We require that La < L^, and thus Ra > Ru'- 
thus there is a value of ha^l,, say hmax, above which the axion flux drops below the neutrino 
flux again. Hence the range hmm < ha^p < hmax is excluded because the axion flux would 
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dominate the one coming from the neutrinos. If one is to remain conservative, many body 
effects should be taken into account, in which case the hmit from this constraint is ||28[ 



ihl + 2hiy/'< 2.85x10-''^. (58) 

This constraint puts bounds on Va, which we write as 

Va > (0.33 X 10^°) {A cos^ 2(3 + B cos 2(3 + Cfl'^ GeV. (59) 

The values of A, B and C, and their uncertainties (which arise from the uncertainties in 
Ag, z, and w) are displayed in Table 1. 

For completeness, we also give the bounds on the DFSZ model originating from the 
electron and photon couplings: 

> 4 X 10^ cos^ /5 GeV, (60) 
Va > 5.1 X 10^ GeV, (61) 

It is clear that the lower bound on the axion scale from the red giant constraint (which 
operates on the axion-electron coupling) is relaxed by a factor of about 2. The axion-photon 
coupling depends on Va/N, as does the mass of the axion and the cosmological constraint in 
inflationary scenarios (in the axion string scenario, = 1/2). These particular constraints 
on the ratio Va/N, are therefore the same for our variant axion models as for DFSZ. 

On the other hand, it is not so straightforward to compare the nucleon constraints. For 
one to see how these bounds affect Va, we plot graphs of Va against P for every model including 
the DFSZ. One can easily check from (Fig. 1) that the lower bound is shifted downwards 
by a factor of 1.4-35 compared to that of the DFSZ axion. We have also displayed how the 
uncertainties in the nucloeon couplings are propagated into uncertainties in the bounds. 



V. CONCLUSIONS 



Although the axion solution to the strong CP problem is one of the most physically 
appealing, axions themselves face a great problem. Despite the fact that they interact very 
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weakly with matter and so are very difficult to track, particle physics experiments together 
with astrophysical considerations and cosmology have managed to constrain axion models 
significantly. In this paper we have found the constraints on for axion models with non- 
standard couplings to quarks and leptons, using data from E143 to determine the values of 
the nucleon couplings, which provide the strongest constraint. We find that the bounds are 
generally weakened, as the nucleon couplings in our variant axion models are smaller. The 
most spectacular effect is for Model I near /3 ~ 7r/4: the bound dips to about Va> 2^ 10^ 
GeV, about a factor 35 less than the DFSZ value. Models I and II have the desirable feature 
that the QCD anomaly coefficient = |, which means that they have no domain wall 
problem, and therefore are viable models for an axion string scenario. For Model I, the 
lower bound on Va can dip to 2 x 10^ GeV for values of (3 near 7r/4. Recalling the upper 
bound on Va in the axion string scenario, Va < 2.6 x 10^^ GeV for Hq = 50 km s~^ Mpc^^ 



1^ , we see that the 'window' for this axion string scenario is actually quite large. 

In an inflationary scenario, the cosmological upper bound is on Va/N, and so the upper 
bound on Va itself is reduced by a factor of 6 for Models I and II, and a factor 3 for Model 
III. 
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TABLES 



Model 


A 


B 


C 


I 


1.13 ±0.14 


0.62 ±0.18 


0.09 ±0.06 


II 


1.13 ±0.14 


2.75 ±0.31 


1.75 ±0.26 


III 


1.13 ±0.14 


1.11 ±0.31 


0.38 ±0.17 


DFSZ 


4.52 ±0.14 


6.42 ± 1.76 


3.75 ± 1.35 



TABLE I. Table of values of coefficients in Eq. (|59|), with experminetal uncertainties, for the 
three variant axion models considered, compared with the values for the canonical DFSZ model. 
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FIGURES 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

P(rad) 

FIG. 1. The lower bound on the Peccei-Quinn symmetry-breaking scale Va for variant axion 
models described in the text and the DFSZ model, plotted as a function of p. Here, tan/3 = V2/V1, 
where vi and V2 are the vacuum expectation values of the two higgs fields giving masses to the 
fermions. The dashed and dot-dashed lines indicate the uncertainties arising from the nucleon 
couplings. 
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